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Typical  microstructures  of  LSM/YSZ,  NiO/YSZ,  Ni/YSZ  composite  electrodes  are  simulated  by  the  dis¬ 
crete  element  method.  The  numerical  microstructures  are  generated  by  taking  into  account  in  a  realistic 
manner  the  sintering  process.  This  allows  complex  microstructures  such  as  bilayers  or  microstructures 
containing  pore  formers  to  be  obtained.  NiO  particles  in  the  NiO/YSZ  composite  electrodes  are  reduced  to 
Ni.  Reduction  is  carried  out  with  the  discrete  element  formalism  which  allows  particle  rearrangement  to 
be  taken  into  account.  We  show  that  the  mechanical  percolation  of  the  YSZ  phase  plays  an  important  role 
during  the  reduction  of  NiO.  The  various  numerical  microstructures  generated  by  sintering  and  reduction 
are  analyzed  to  evaluate  important  microstructural  features  such  as  macroscopic  porosity,  pore  surface 
area  and  Triple  Phase  Boundary  length. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electrochemical  performance  of  an  SOFC  (Solid  Oxide  Fuel  Cell) 
or  SOEC  (Solid  Oxide  Electrolysis  Cell)  strongly  depends  on  the 
microstructure  of  the  porous  electrodes  which  constitute  the  most 
critical  element  of  the  cell.  Electrodes  are  fabricated  from  ceramic 
powders  that  are  partially  sintered  to  obtain  a  functionalized 
porous  ceramic.  Taking  advantage  of  the  freedom  in  microstructure 
design  that  powder  processes  allow,  experimentalists  have  pro¬ 
posed  various  microstructural  architectures  for  porous  electrodes. 
They  have  well  established  the  strong  dependence  of  the  functional 
properties  of  a  porous  ceramic,  at  a  fixed  chemical  composition,  on 
microstructure  [1-6]. 

Thus,  a  number  of  studies  have  demonstrated  that  gains  may  be 
obtained  by  improving  the  microstructure  of  the  porous  electrodes 
for  a  given  cathode  or  anode  material.  The  main  goal  of  improving 
the  microstructural  design  is  to  lower  the  operating  temperature 
without  decreasing  the  power  densities.  For  example,  the  fuel  cell 
program  of  Forschungszentrum  Jiilich  developed  so-called  func¬ 
tional  layers,  both  at  the  cathode/electrolyte  interface  and  at  the 
anode  substrate/electrolyte  interface  [3].  In  that  case  the  cathode 
functional  layer  is  a  composite  made  of  a  mixture  of  electronic 
and  ionic  conducting  particles  which  allows  a  large  increase  of  the 
number  of  active  sites.  The  use  of  fine-grained  microstructure  also 


*  Corresponding  author.  Tel.:  +33  476826337;  fax:  +33  476826382. 
E-mail  address:  Christophe.Martin@grenoble-inp.fr  (C.L.  Martin). 

0378-7753 /$  -  see  front  matter  ©  2010  Elsevier  B.V.  All  rights  reserved, 
doi:  10.1 01 6/j.jpowsour.201 0.09.033 


increases  the  number  of  so-called  Three-Phase  Boundary  (TPB)  per 
unit  volume.  Typically  submicron  sized  powders  (provided  they 
do  not  experience  too  large  a  coarsening  during  sintering)  are 
expected  to  produce  larger  TPB  length  per  unit  volume  of  electrode 
[7,8].  Other  examples  of  enhanced  microstructures  deal  with  the 
control  of  the  pore  architecture  both  for  anodes  and  cathodes  [5,9] 
and  with  the  possibility  of  grading  porosity  along  the  thickness  of 
the  electrodes  [  1 0,1 1  ].  Control  of  porosity  and  pore  architecture  has 
been  demonstrated  with  various  types  of  pore  formers  in  electrode 
fabrication  [5,8,9,12].  Pore  formers  are  used  to  adjust  the  shrink¬ 
age  of  each  layer  during  the  co-sintering  stage  of  multi-layered 
electrodes  but  may  also  be  used  to  enhance,  under  operation,  the 
diffusion  of  gaseous  species  to  the  functional  layers  [12,13]. 

These  examples,  using  very  different  routes  for  optimizing 
electrodes,  show  that  it  is  critical  to  gain  a  more  thorough  under¬ 
standing  of  the  correlation  between  the  microstructure  at  the 
length  scale  of  powder  particles  and  the  macroscopic  electrode  per¬ 
formances.  Also,  linking  pertinent  microstructural  features  to  the 
macroscopic  behavior  would  allow  predictive  models  to  guide  the 
design  of  optimized  microstructural  architectures. 

Since  most  routes  for  fabricating  electrodes  involve  deposit¬ 
ing  a  powder  mixture  onto  a  dense  electrolyte  and  then  sintering 
this  mixture,  such  models  should  be  able  to  reproduce  these  steps 
with  some  accuracy.  This  would  allow  mimicking  the  complex 
microstructural  designs  produced  by  powder  processing.  How¬ 
ever,  such  a  modeling  approach  requires  an  integrated  approach 
of  the  electrode  processing  from  the  initial  powders  to  the  final 
porous  electrode.  Thus,  any  model  of  the  electrode  microstructure 
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needs  to  acknowledge  the  particulate  nature  of  the  material  which 
has  been  processed  to  obtain  the  final  electrode.  This  has  been 
only  partially  recognized  by  analytical  models  which  compute  for 
example  TPB  length  as  a  function  of  porosity,  particle  size  and  com¬ 
position  [14-17].  These  analytical  models  have  the  advantage  of 
demonstrating  in  a  simple  and  elegant  way  the  effect  of  the  main 
parameters  which  control  the  TPB  length.  However,  the  microstruc¬ 
ture  is  still  modeled  in  a  very  approximate  manner  and  the  sintering 
process  is  treated  only  through  its  effect  on  the  decrease  of  the 
macroscopic  porosity. 

Other  models  use  numerical  microstructures  to  mimic  the 
porous  electrodes.  These  numerical  simulations  use  a  random  pack¬ 
ing  of  spheres  which  indent  each  other  to  create  a  porous  material. 
The  contacts  between  electronic  and  ionic  conducting  particles 
allow  the  TPB  length  to  be  computed  for  various  compositions 
and  electrode  thicknesses  [18,19],  particle  sizes  [20]  and  particle 
size  distributions  [21  ]  and  composition  gradients  [22].  Using  image 
analysis  techniques,  it  is  also  possible  to  obtain  from  these  numeri¬ 
cal  microstructures  other  important  microstructural  features  such 
as  pore  size  and  geometry  and  pore  percolation  threshold  [21,23]. 
From  these  microstructural  features  and  from  the  intrinsic  conduc¬ 
tive  properties  of  the  particles,  it  is  then  possible  to  compute  the 
effective  electric  and  electrochemical  properties  of  the  electrodes 
by  assuming  that  current  is  conducted  at  solid  bonds  between  par¬ 
ticles  [14,17,18,20,22,24,25]. 

All  these  models  recognize  the  particulate  (or  discrete)  nature  of 
the  porous  electrode  and  have  improved  our  understanding  of  the 
effect  of  important  features  such  as  percolation,  electrode  thick¬ 
ness,  composition  and  particle  size.  However,  they  are  still  a  very 
approximate  picture  of  the  real  microstructure.  In  particular  in 
most  of  these  models,  the  sintering  process  is  treated  simply  as 
the  growth  of  particle  radii  to  form  solid  bonds  between  parti¬ 
cles  [17,20,21,24,25],  thus  neglecting  particle  rearrangement  that 
must  occur  at  interfaces  with  the  dense  electrolyte  or  when  a  mix¬ 
ture  of  particles  with  different  sizes  is  used.  Also,  the  reduction 
of  the  NiO  ceramic  into  metallic  Ni  to  form  a  cermet  anode  is  not 
treated  in  terms  of  its  microstructural  consequences.  For  example 
Golbert  et  al.  have  treated  the  reduction  of  NiO  and  the  result¬ 
ing  increase  in  porosity  by  introducing  numerically  pore  formers 
which  are  subsequently  removed  [20].  The  most  advanced  models 
on  reduction  of  NiO  to  Ni  (and  subsequent  reoxidations)  simplify 
the  microstructure  to  a  continuum,  where  stresses  and  strains  are 
one-dimensional  [26]. 

In  that  context,  the  aim  of  this  work  is  to  provide  more  realistic 
numerical  microstructures  that  better  reflect  electrode  processing. 
In  particular,  the  sintering  process  should  be  explicitly  taken  into 
account  to  better  reflect  the  complexity  of  the  resulting  microstruc¬ 
ture.  Using  the  discrete  element  method,  which  enforces  force 
equilibrium  for  each  particle  in  the  packing  during  sintering,  we 
show  that  indeed  it  is  possible  to  reproduce  complex  microstruc¬ 
tures  with  pore  formers  and/or  multi-layers. 

In  typical  SOFC/SOEC’s,  02  electrodes  are  bilayers  usually  made 
of  yttria-stabilized  zirconia  (YSZ)  and  perovskite-type  oxide  like 
Lai_xSrxMn03  (LSM).  For  state-of-the-art  H2  electrodes,  compos¬ 
ites  like  Ni/YSZ  or  Ni/cerium-gadolinium  oxide  (CGO)  cermets  are 
popular  candidates.  The  simulation  method  is  thus  applied  here  to 
typical  LSM/YSZ  and  NiO/YSZ  (before  reduction  to  Ni)  composite 
electrodes.  In  any  case,  the  simulation  method  is  general  enough 
to  be  applied  to  other  materials.  For  example,  although  CGO  par¬ 
ticles  have  mixed  conducting  properties,  they  may  be  handled  by 
the  same  method  since  they  should  sinter  similarly  to  YSZ  and  they 
have  very  similar  elastic  properties  [27]. 

The  paper  is  organized  as  follows:  In  Section  2,  we  present 
the  particle  packings  that  are  used  for  generating  electrodes  with 
pore  formers  and/or  layers.  Sintering  of  the  powder  packings  is 
described  in  Section  3.  The  reduction  of  NiO  into  Ni  is  simulated  in 


Fig.  1.  Section  of  a  NiO/YSZ-pf  (with  pore  formers)  electrode;  (a)  before  sintering 
and;  (b)  after  sintering.  Pore  formers  are  removed  before  sintering.  Colors  indicate 
the  number  of  contacts  Z  for  each  particle.  (For  interpretation  of  the  references  to 
color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 

Section  4  while  the  microstructural  characteristics  of  the  numerical 
microstructures  are  described  and  commented  in  Section  5.  When¬ 
ever  possible,  we  attempt  to  compare  our  results  with  experimental 
data  from  the  literature. 

2.  Porous  electrodes  generation 

Typical  H2  and  02  electrodes  used  in  the  SOEC/SOFC  technolo¬ 
gies  [5]  have  been  considered  in  this  work.  It  should  be  noted 
that  porous  electrodes  can  be  fabricated  by  different  routes  (tape 
casting,  screen  printing,  . . .)  depending  on  the  experience  of  the 
manufacturers  and  the  final  purpose  of  the  electrodes.  However, 
it  is  out  of  the  scope  of  this  paper  to  detail  each  step  of  such  pro¬ 
cesses.  Instead,  we  only  considered,  for  the  sake  of  simplicity,  that 
electrodes  are  deposited  in  a  green  state  and  then  consolidated 
by  sintering  to  obtain  their  final  microstructure.  According  to  our 
own  experience,  typical  green  electrodes  are  basically  constituted 
by  a  loose  packing  of  active  materials  particles  (NiO,  YSZ,  LSM) 
and  of  pore  formers  whom  characteristics  (size  and  volume  frac¬ 
tion)  are  summarized  in  Table  1 .  Any  additional  precursors  of  these 
microstructures  are  not  taken  into  account  in  this  study. 

Rectangular  green  electrodes  were  generated  numerically  by 
using  the  procedure  summarized  in  Fig.  1  and  described  hereafter. 
A  packing  of  spherical  particles  located  at  random  positions  is  first 
generated  with  the  constraint  that  particles  should  not  contact  each 
other.  This  initial  packing  is  in  a  simulation  box  bounded  by  rigid 
planes  in  the  z  direction  and  by  periodic  conditions  in  the  x  and  y 
directions.  At  this  initial  stage,  particle  size  distribution  and  vol¬ 
ume  fractions  of  each  type  of  particle  are  defined  depending  on  the 
final  purpose  of  the  electrode.  The  particle  packing  is  then  densified 
with  stiff  elastic  interactions  between  particles  to  obtain  a  dense 
green  body.  The  procedure  to  generate  such  packings  of  particles 
has  been  explained  in  previous  works  [28-30].  Fig.  la  shows  the 
green  microstructure  at  this  stage. 

The  first  type  of  H2-electrode,  denoted  as  NiO/YSZ  is  obtained 
from  a  mixture  of  NiO  and  YSZ  particles  of  the  same  size 
(0.5  pan  ±5%).  By  default,  the  volume  percentage  of  NiO  in  the  mix¬ 
ture  is  65%.  The  second  type  of  H2 -electrode,  made  of  the  same 
NiO/YSZ  mixture  as  for  the  first  electrode,  contains  additionally 
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Table  1 

Electrode  characteristics,  pf  indicates  electrodes  with  pore  formers.  The  LSM/YSZ-pf  electrode  is  a  bilayer  made  of  a  YSZ  +  LSM  layer  and  of  a  LSM  layer.  Two  volume  percentage 
of  NiO  have  been  studied  for  the  NiO/YSZ  and  NiO/YSZ-pf  electrodes. 


NiO/YSZ 

NiO/YSZ-pf 

LSM/YSZ-pf 

Particle  size  (±5%) 

NiO^O.5  pan 

NiO^O.5  pan 

LSM^l  pm 

YSZ^0.5  pm 

YSZ^0.5  pm 

YSZ^0.5  pm 

Volume  fraction 

65%  NiO  and  35%  NiO 

65%  NiO  and  35%  NiO 

YSZ  +  LSM— >-50%  YSZ 

LSM^0%  YSZ 

Pore  former  volume  fraction 

32% 

YSZ  +  LSM^-32% 

LSM^30% 

Pore  former  size  (±5%) 

2  (Jim 

YSZ  +  LSM— >-1.4  pm 
LSM^2.4  pm 

Sintered  electrode  height 

30  pm 

40  pan 

YSZ  +  LSM^30pm 
LSM— >-23  pm 

Number  of  particles 

10,000 

73,000 

92,000 

fugitive  pore  formers  (2  p,m  ±5%).  It  is  denoted  as  NiO/YSZ-pf.  Once 
pore  formers  are  removed  from  this  second  electrode,  the  volume 
percentage  is  65%  volume  of  NiO  as  for  the  first  electrode.  Note  that 
the  electrode  is  not  graded  and  that  the  active  zone  is  assumed 
to  spread  all  along  the  thickness.  The  microstructure  of  this  H2- 
electrode  is  illustrated  in  Fig.  lb  before  the  reduction  step.  The 
partially  sintered  microstructure  exhibits  two  populations  of  pores: 
the  small  pores  left  in  between  partially  sintered  particles,  and  the 
large  pores  originating  from  the  removal  of  the  pore  formers.  This 
double  porosity  compares  qualitatively  well  with  the  anodes  sin¬ 
tered  with  fugitive  pore  formers  by  Suzuki  et  al.  [5].  These  anodes 
were  made  of  NiO-Sc-stabilized  zirconia  and  Ce-doped  zirconia. 
The  double  porosity  is  linked  to  the  relatively  low  temperature  used 
by  these  authors  (1250  °C),  which  ensures  only  partial  sintering  of 
the  particles,  and  is  clearly  observed  in  their  Fig.  2  for  cell  A. 

The  third  electrode  is  a  02  electrode  made  of  two  layers  which 
are  subsequently  co-sintered  together.  The  active  layer,  adjacent 
to  the  electrolyte,  is  a  mixture  of  YSZ  (0.5  fxm±5%)  and  of  LSM 
(l|jim±5%)  particles  and  of  pore  formers  (1.4  pum  ±5%),  while 
the  diffusion  layer  is  made  of  LSM  particles  (1  pan  ±5%)  with 
2.4  pm  ±5%  fugitive  pore  formers.  This  electrode  is  denoted  as 
LSM/YSZ-pf.  These  electrodes  are  typical  of  composite  electrodes 
encountered  in  the  literature  [5].  A  typical  microstructure  of  such  a 
bilayered  electrode  observed  after  the  sintering  stage  on  commer¬ 
cial  cells  is  illustrated  in  Fig.  2.  The  narrow  size  distribution  that 
we  have  used  is  obviously  an  oversimplification  of  real  microstruc¬ 
tures.  Literature  provides  examples  of  porous  electrodes  with  a 
rather  uniform  size  distribution  [3-5,7]  (and  also  Fig.  2  in  this  work) 
and  examples  where  the  size  distribution  is  much  larger  [8,9,13]. 
For  sake  of  simplicity,  we  have  chosen  to  leave  the  issue  of  large 
particle  size  distribution  for  future  work. 


LSM 


YSZ+ 

LSM 


Fig.  2.  Microstructure  of  the  LSM/YSZ  bilayer.  The  white  dotted  curve  indicates  the 
interface  between  the  composite  layer  YSZ  +  LSM  and  the  homogeneous  LSM  layer. 
YSZ  and  LSM  cannot  be  distinguished  in  this  SEM  micrograph. 


Fig.  la  shows  the  microstructure  of  the  NiO/YSZ-pf  pack¬ 
ing  before  sintering.  The  large  pore  former  particles  are  packed 
together  with  the  other  particles.  Colors  indicate  the  coordination 
number  (number  of  contacts  with  other  particles)  for  each  indi¬ 
vidual  particle.  The  coordination  number  is  a  good  indicator  of  the 
local  packing  density.  For  electrodes  containing  pore  formers  such 
as  the  electrode  depicted  in  Fig.  1,  pore  former  particles  are  elimi¬ 
nated  before  the  numerical  sintering.  This  is  to  mimic  real  ceramic 
processes  in  which  fugitive  polymer  pore  formers  are  introduced 
and  then  burned  out  at  high  temperature  [5,8,9].  All  packings  have 
£0  =  0.5  initial  porosity  after  fugitive  pore  former  removal,  prior  to 
sintering. 

Electrode  heights  are  typical  of  electrolyte-supported  cells  (30 
to  40  |jim).  Also,  it  should  be  noted  that  in  order  to  keep  reasonable 
CPU  time  for  the  simulation,  it  is  not  possible  to  consider  the  whole 
width  of  the  electrode  in  the  x  and  y  directions.  Instead,  we  simu¬ 
late  a  pillar  which  is  typically  15  to  20  pm  wide.  Still,  the  number 
of  simulated  particles  (10,000  for  the  NiO/YSZ  electrode  without 
pore  former  and  >70,000  for  the  electrodes  with  pore  formers) 
should  be  large  enough  to  ensure  that  the  numerical  microstructure 
constitutes  a  representative  volume  element  of  the  electrode. 


3.  Sintering 


All  packings  are  subjected  to  a  numerical  sintering  procedure 
which  has  been  detailed  elsewhere  [30-32].  Briefly,  it  consists  in 
introducing  contact  forces  that  account  for  sintering  phenomena 
between  particles  at  high  temperature.  Bouvard  and  McMeeking’s 
model,  which  considers  grain  boundary  and  surface  diffusion  to  be 
the  main  mechanisms  of  mass  transport,  is  used  [33].  We  denote 
ys  as  the  surface  energy  and  Ab  as  the  diffusion  parameter: 

^b=kf^bDbl  (1) 

where  Db  =  Dobexp  ( -  QblRT)  is  the  diffusion  coefficient  for  vacancy 
transport  in  the  grain  boundary  with  thickness  8b  and  activation 
energy  and  £2  is  the  atomic  volume.  For  two  spherical  particles 
of  radii  R i  and  R2 ,  having  a  contact  radius  ab  and  indentation  hb 
(Fig.  3),  the  normal  force  Ns  acting  on  the  contact  is  given  by: 


no*  dht 
5  8Ab  dt 


9 

4 


7TR*Ys, 


(2) 


where  ft*  =  ftift2/(l?i  +fi2).  The  normal  contact  force  in  Eq.  (2)  is 
the  sum  of  two  terms.  The  first  term  is  a  normal  viscosity  term, 
which  acts  both  in  tension  or  compression  depending  on  the  rate 
of  approach  of  the  two  particles  (dhbldt).  The  second  term,  linked 
to  surface  energy,  acts  only  in  tension  and  is  independent  of  contact 
size. 

Similarly,  a  tangential  contact  force  Ts  opposes  the  tangential 
component  of  the  relative  velocity  at  the  contact,  du/dt,  and  is  given 
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Fig.  3.  Schematic  of  two  particles  forming  a  sintered  bond  of  radius  Qj,.  The  Scapi 
surface  is  subtracted  to  the  particle  surface  while  Scy/  is  added  to  compute  the  total 
free  surface  in  Eq.  (6). 


by  [34,35]: 


na2R*2 

"lATdt’ 


(3) 


where  r]  is  a  viscous  parameter  with  no  dimension.  In  the  simula¬ 
tions,  ij  =  10-3  for  contacts  between  particles  and  r]  =  0  for  contacts 
between  particles  and  the  substrates  [32].  We  assume  that  con¬ 
tacts  become  rapidly  large  enough  to  oppose  any  rotation  of  the 
particles. 

For  sintering,  free  surfaces  along  the  x  and  y  directions  are  used 
and  replace  the  periodic  conditions  previously  used  for  compaction. 
The  two  rigid  planes  with  normal  along  the  z  direction  are  main¬ 
tained  during  sintering  and  thus  act  as  planar  substrates.  Particles 
form  sintered  bonds  both  with  neighboring  particles  and  with  the 
two  substrates.  The  bond  size  ab  between  two  particles  of  radii  R i 
and  R2  grows  according  to  Coble’s  model  [36]: 


d ab  _  2  R-  dhb 
d t  ~  ab  dt  ’ 


which  is  in  good  accordance  with  various  numerical  simulations 
on  pairs  of  particles  of  equal  size  [33]  or  different  sizes  [37].  When 
integrating  Eq.  (4)  from  an  initial  value  ab  =  0,  it  yields  a  bond  radius 
a 2  =  4 R*hb  which  is  larger  (by  a  V2  factor)  than  the  contact  radius 
obtained  from  a  simple  geometric  intersect  of  two  overlapping 
spheres  (ac).  This  is  because  Coble’s  model  accounts  for  volume  con¬ 
servation  as  sintering  proceeds.  This  larger  contact  size  impacts  on 
the  calculation  of  the  Triple  Phase  Boundary  length  and  is  neglected 
by  models  that  approximate  sintering  to  the  growth  of  particle  radii 
[17,20,21,24,25]. 

As  shown  by  Fig.  3,  sintering  results  in  particle  centers  approach¬ 
ing  each  other  and  thus  particles  becoming  truncated  spheres.  In 
other  words,  densification  of  the  powder  proceeds  in  the  discrete 
element  simulations  by  the  gradual  overlap  of  particles.  This  is 
accompanied  also  by  the  appearance  of  new  contacts  between  par¬ 
ticles.  This  is  clearly  demonstrated  by  comparing  Fig.  1  a  and  Fig.  1  b, 
which  also  shows  a  decrease  of  the  height  of  the  particulate  assem¬ 
bly. 

There  is  a  lack  of  data  corresponding  to  the  sintering  material 
parameters  relevant  for  SOFC  manufacturing.  Thus,  we  do  not  have 
access  to  the  sintering  material  parameters  that  are  needed  in  Eqs 
(1 ),  (2)  and  (3)  for  each  type  of  ceramic  particle  modeled  here  (NiO, 
YSZ  and  LSM).  To  tackle  this  problem,  we  make  the  assumption  that 
all  particles  of  a  given  size  sinter  with  the  same  kinetics,  irrespec¬ 
tive  of  their  composition.  Under  this  assumption  and  recognizing 
that  the  exact  kinetics  of  densification  of  the  porous  microstructure 


(the  time  necessary  to  attain  a  given  relative  density  at  a  given  tem¬ 
perature)  is  out  of  reach,  we  may  obtain  realistic  microstructures. 

This  is  because  examination  of  Eq.  (2)  indicates  that  for 
a  r  =  R4l(Abys),  where  R  is  a  length  scale  equal  to  the  mean 
particle  radius,  gives  a  normalized  time  which  introduces  time- 
temperature  equivalence  (through  the  Ab  term).  The  expression 
of  r  is  obtained  by  enforcing  zero  contact  force  (force  equilibrium, 
Ns  =  0)  and  by  relating  lengths  ( R *,  ab  and  hb )  to  the  mean  radius 
R.  This  normalized  time  correctly  reveals  that  large  particles  sin¬ 
ter  slower,  and  that  larger  temperature  are  associated  with  shorter 
sintering  time.  Using  this  normalized  time  in  our  simulation,  we 
may  use  an  arbitrary  set  of  material  parameters.  Here,  we  have 
used  alumina  parameters  ( Dob  =  1 .3 1 O-08  m3  s-1 ,  £2  =  8.47 1 O-30  m3 , 
Qb  =  475  kj  mole-1 ,  and  ys  =  1 .1 J  m-2 )  [38,39].  Thus,  while  the  phys¬ 
ical  time  necessary  to  reach  a  given  microstructure  is  dependent  on 
this  set  of  material  parameters  (and  on  the  sintering  temperature), 
the  microstructure  itself  is  independent  of  the  material  parameters 
in  our  model. 

However,  differences  in  particle  size  in  the  packing  play  an 
important  role  in  the  sintering  model  through  the  R*  effective 
radius.  In  particular,  Eqs.  (2)  and  (4)  show  that  the  time  necessary 
to  decrease  porosity  to  a  given  value  in  a  volume  element  made  of 
particles  of  size  R  scales  with  R-4  [37,40].  Thus,  when  force  equi¬ 
librium  is  ensured,  the  rate  of  approach  (dh^/dt)  between  two  YSZ 
particles  (0.5  pm  size)  or  between  a  YSZ  particle  and  an  LSM  particle 
(1  p,m  size)  is  larger  than  between  two  LSM  particles. 

Note  also  that  for  simplicity  we  have  not  considered  coarsening 
between  particles  during  sintering  as  in  a  previous  work  [31  ].  This 
is  a  reasonable  assumption  for  the  partial  sintering  process  which 
we  are  modeling  here. 

The  sintering  procedure  presented  above  has  been  applied  to  the 
three  electrodes  described  in  Table  1 .  Starting  from  an  initial  poros¬ 
ity  of  £0  =  0.5,  sintering  was  interrupted  at  decreasing  amounts  of 
porosity  to  study  the  evolution  of  the  main  microstructural  features 
which  are  pertinent  to  electrochemical  processes  (see  Section  5). 

Concerning  the  cosintering  of  the  bilayered  LSM/YSZ-pf  elec¬ 
trode,  Fig.  4  shows  the  evolution  of  the  local  porosity  s  along  the 
z  axis  at  various  sintering  times.  The  two  layers  have  started  from 
the  same  approximate  initial  porosity  (s0  ~  0.5).  The  initial  porosity 
exhibits  some  deviation  from  the  0.5  mean  value  due  to  the  finite 
size  of  the  simulated  volume.  Also,  the  initial  interface  between  the 
LSM  layer  and  the  LSM  +  YSZ  layer  exhibits  a  smaller  porosity.  This  is 
due  to  the  fact  that  by  construction  of  the  bilayer,  no  large  pore  can 
cross  the  interface.  Recall  that  large  pores  are  initially  introduced 
in  the  packing  as  pore  former  particles  (see  Fig.  1 ). 

Note  how  this  denser  initial  interface  has  translated  from  32  pm 
height  to  27  p,m  height  and  has  kept  its  relatively  lower  porosity  as 
sintering  proceeded.  More  importantly,  because  the  YSZ  particles 
are  smaller  (0.5  pan)  than  the  LSM  particles  (1  pm),  the  YSZ-LSM 
layer  sinters  faster  than  the  LSM  layer.  The  larger  densification  of 
the  YSZ-LSM  layer  is  observable  from  the  surface  distortion  at  the 
interface  with  the  LSM  layer  that  constrains  somewhat  the  bottom 
layer.  Recall  that  the  bilayer  exhibits  free  surfaces  in  the  x  and  y 
directions  during  the  sintering  stage.  At  a  mesoscopic  length  scale, 
curve  b  in  Fig.  4  clearly  indicates  that  the  porosity  of  the  YSZ-LSM 
layer  becomes  much  smaller  ( s  ^  0.1 8)  than  the  porosity  of  the  LSM 
layer  ( s  ^  0.35)  as  sintering  proceeds. 

This  is  an  encouraging  indication  that  our  sintering  model 
is  behaving  correctly.  Indeed,  we  have  observed  on  experimen¬ 
tally  cosintered  bilayered  LSM/YSZ  electrodes  (Fig.  2)  with  initial 
microstructures  characterized  by  the  properties  given  in  Table  1 
that  the  average  porosity  of  the  LSM  and  of  the  YSZ-LSM  layers 
are  0.34  and  0.18,  respectively.  These  numbers  compare  favorably 
with  curve  b  in  Fig.  4.  Again,  it  should  be  clear  that  the  exact  heating 
cycle  (time  and  temperature)  necessary  to  obtain  the  microstruc¬ 
ture  defined  by  curve  b  cannot  be  reproduced  by  our  model  because 
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Fig.  4.  Evolution  of  the  porosity  along  the  z  axis  of  the  LSM/YSZ-pf  electrode  (a)  before  sintering  and  (b)  after  sintering.  The  corresponding  final  bilayered  microstructure  is 
shown  on  the  right  hand  side.  Arrows  indicate  the  experimentally  measured  average  porosity  for  each  layer. 


it  would  necessitate  the  knowledge  of  diffusion  coefficients  and 
surface  energies  which  are  not  available  for  YSZ  and  LSM.  Still,  if  the 
sintering  simulation  is  stopped  when  the  porosity  of  the  YSZ-LSM 
layer  is  approximately  0.18,  the  model  appropriately  predicts  that 
the  LSM  layer  porosity,  which  contains  larger  particles,  is  approxi¬ 
mately  0.35. 


4.  Reduction  of  NiO  Ni 

Once  sintered,  the  NiO  ceramic  particles  in  the  anode  are 
reduced  to  Ni  metal  on  exposure  to  hydrogen.  The  reduction  of  the 
NiO  particles  has  been  numerically  simulated  by  a  simple  proce¬ 
dure  based  on  the  model  proposed  by  Klemenso  et  al.  [41  ].  The  YSZ 
particles  are  bonded  together  with  stiff  elastic  bonds  (Young’s  mod¬ 
ulus  220  GPa)  and  infinite  strength.  The  NiO  particles  are  bonded 
together  and  to  YSZ  particles  by  weaker  bonds  that  may  deform 
easily  if  strained  (Young’s  modulus  2  GPa).  This  is  to  model  the  fact 
that  the  NiO  ceramic  particles  become  softer  metallic  particles  at 
high  temperature  during  reduction.  Still,  it  should  be  clear  that  we 
do  not  attempt  to  model  the  plastic  (or  viscoplastic)  deformation 
that  Ni  particles  may  undergo  at  high  temperature.  The  NiO  parti¬ 
cles  are  then  submitted  to  a  gradual  diameter  reduction  to  account 
for  the  NiO  Ni  phase  change.  We  have  chosen  to  look  at  two 
values  for  the  volume  reduction  of  the  nickel  phase:  25%  and  41%, 
corresponding  to  9%  and  16%  radius  reduction,  respectively.  The 
value  of  25%  may  be  regarded  as  typical  of  the  amount  of  a  partial 
reduction  [41-43],  while  the  41%  is  the  upper  theoretical  limit  for 
the  full  reduction  of  NiO  to  Ni  [44].  The  size  reduction  implies  a 
reorganization  of  the  Ni  and  YSZ  particles  which  is  accounted  for 
in  the  discrete  simulations  by  imposing  force  equilibrium  to  each 
particle.  Note  that,  contrarily  to  the  sintering  stage,  the  sample  is 
not  bound  to  any  substrate  in  the  reduction  stage. 

Two  NiO  volume  fractions  (ONi0  =  35%  and  ONi0  =  65%)  were 
simulated  to  study  the  effect  of  mechanical  percolation  of  the  YSZ 
phase.  By  mechanical  percolation  we  mean  that  a  rigid  skeleton 
of  YSZ  particles  that  can  sustain  forces  without  large  deformation 
has  formed.  Such  percolation  is  always  larger  than  the  more  classi¬ 
cal  geometrical  percolation  threshold  [45].  When  the  YSZ  volume 
fraction  is  only  35%  (ONi0  =  65%),  mechanical  percolation  is  not 
enforced.  This  has  some  important  effect  on  the  porosity  increase 
that  should  be  brought  by  NiO  reduction.  On  the  contrary,  when  the 
YSZ  volume  fraction  is  65%  (ONi0  =  35%),  mechanical  percolation  is 
ensured.  In  that  case,  owing  to  the  large  rigidity  of  the  YSZ-YSZ 
bonds,  a  YSZ  skeleton  is  formed  and  it  is  mostly  the  Ni  particles 
that  rearrange  during  reduction. 

If  the  YSZ  phase  percolates  mechanically,  it  may  be  assumed 
that  no  macroscopic  shrinkage  will  occur  upon  reduction  since  a 


YSZ  skeleton  holds  the  sample.  In  that  case,  the  porosity  s  may  be 
simply  deduced  from  the  following  equation: 

8  =  £i  +  AI/Nio^Ni(l  -  8i)<& NiO>  (5) 

where  £z-  is  the  porosity  before  reduction  and  AVNio^Ni  is  the  vol¬ 
ume  reduction  of  the  Ni  phase  (from  0  to  0.41 ).  Fig.  5a  and  b  shows 
the  evolution  of  porosity  as  reduction  proceeds  (AVNio^Ni  =  0  to 
41%)  for  samples  with  and  without  pore  formers,  respectively. 
They  indicate  that  indeed,  when  mechanical  percolation  is  enforced 
(ONi0  =  35%),  Eq.  (5)  is  in  good  accordance  with  simulations.  This  is 
consistent  with  experimental  data  that  report  very  little  macro¬ 
scopic  shrinkage  upon  the  first  reduction  for  materials  which 
contain  a  sufficient  amount  of  YSZ  phase  [41,46,47]. 

However,  when  the  NiO  volume  fraction  is  high  (ONi0  =  65%), 
the  porosity  increase  given  by  Eq.  (5)  overestimates  the  porosity 
increase  calculated  from  the  simulations.  This  is  because,  in  that 
case,  the  sample  experience  macroscopic  shrinkage.  Experimen¬ 
tal  data  confirms  that  in  certain  cases,  macroscopic  shrinkage  may 
occur,  thus  limiting  the  porosity  increase  of  the  sample  [46,48].  It 
has  also  been  reported  that  significant  shrinkage  may  occur  when 
the  volume  fraction  of  YSZ  is  small  (ONi0  =  90%)  [48].  This  is  in  good 
accordance  with  the  present  simulations  which  demonstrate  the 
importance  of  the  formation  of  a  YSZ  skeleton  upon  reduction. 

The  curves  shown  in  Fig.  5a  and  b  are  reported  for  an  initial 
porosity  {sit  porosity  before  reduction)  of  approximately  30%.  The 
same  qualitative  behavior  was  observed  for  the  other  initial  porosi¬ 
ties. 

In  accordance  with  the  porosity  evolution  reported  above,  we 
have  observed  different  reduced  microstructures,  depending  on 
the  initial  YSZ  volume  fraction  (35%  or  65%).  If  the  YSZ  volume 
fraction  is  large  (65%),  the  resulting  microstructure  accommodates 
the  volume  change  predominantly  by  pore  opening  inside  the  YSZ 
skeleton  while  macroscopic  volume  change  does  not  occur.  If  the 
volume  fraction  is  small  (35%),  shrinkage  partially  accommodates 
the  volume  change  of  Ni  particles. 

Note  that  by  the  end  of  the  simulated  reduction  procedure 
the  cermet  is  in  the  immediate  reduced  state  as  opposed  to  the 
long-term  reduced  state  which  would  involve  coarsening  of  the 
Ni  particles  [41  ].  Ni  coarsening  is  not  included  in  the  model  in  its 
present  state. 

5.  Microstructural  characterization 

Once  sintered,  the  numerical  microstructures  may  be  char¬ 
acterized  to  obtain  some  useful  information  on  microstructural 
parameters  that  are  critical  for  electrochemical  efficiency.  We  have 
focused  on  two  parameters:  the  pore  surface  area  which  affects 
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Fig.  5.  Evolution  of  the  macroscopic  porosity  of  samples  without  pore  former  (a)  and  with  pore  formers  (b)  upon  reduction  for  two  initial  volume  fractions  of  NiO  (ONi0  =35 
and  65%).  Broken  lines  indicate  the  evolution  of  porosity  when  no  macroscopic  shrinkage  occurs  (Eq.  (5)). 
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Fig.  6.  Evolution  of  the  pore  surface  per  unit  volume  versus  the  macroscopic  poros¬ 
ity  e  of  NiO-YSZ  electrodes  with  and  without  pore  former  before  reduction  and 
of  bilayered  electrode  LSM/YSZ-pf.  The  arrowed  lines  connect  the  LSM  and  the 
YSZ  +  LSM  layers  for  the  initial  (a)  and  final  (b)  microstructures  as  shown  in  Fig.  4. 


interestingly,  Fig.  6  indicates  that  for  a  given  macroscopic  poros¬ 
ity,  introducing  pore  formers  leads  to  a  decrease  of  the  total  free 
surface.  This  is  because  there  are  two  populations  of  pores  induced 
by  the  microstructure  with  pore  formers:  the  natural  small  pores 
left  in  between  particles,  and  the  large  pores  that  are  left  after  the 
removal  of  fugitive  pore  formers.  When  considering  the  total  pore 
surface  per  unit  volume  for  a  given  macroscopic  porosity,  the  sim¬ 
ulations  indicate  that  it  is  more  beneficial  to  use  small  pores  than 
large  ones.  However,  large  percolating  pores  will  certainly  have  an 
overall  beneficial  effect  when  electrode  performance  is  limited  by 
gas  diffusion. 

The  same  analysis  is  carried  out  on  the  LSM/YSZ  +  LSM  bilayer.  In 
that  case,  we  distinguish  the  two  layers  LSM  and  YSZ  +  LSM  in  Fig.  6. 
Because  the  LSM  particles  are  double  the  size  of  the  NiO  or  YSZ  par¬ 
ticles,  the  free  surfaces  per  unit  volume  of  the  LSM  and  YSZ  +  LSM 
layers  are  smaller.  More  generally,  it  may  be  shown  that  for  a  given 
microstructure,  S  is  approximately  inversely  proportional  to  the 
average  particle  size.  Thus,  the  results  of  Fig.  6  may  easily  be  gener¬ 
alized  to  other  particle  sizes.  Also,  for  a  monomodal  powder  made 
of  particles  of  size  2 R,  S  is  bounded  by: 


adsorption  and  surface  diffusion,  and  the  total  length  ofTriple  Phase 
Boundaries  (Ltpb)  where  electrolyte,  electrocatalyst  and  gas  phase 
react.  The  reaction  rate  of  the  three  phases  is  believed  to  be  corre¬ 
lated  to  the  TPB  length. 

5.1.  Pore  surface  area  per  unit  volume  of  sample 

The  pore  surface  area  (free  surface)  per  unit  volume,  is  calcu¬ 
lated  from  the  partially  sintered  assembly  of  N  spherical  particles 
by  considering  their  surface  and  the  surface  geometry  of  each  of 
the  I<  contacts  in  the  electrode  (Fig.  3): 

(N  I< 

E4jrRi-E(s*  ,capl  + 

_  /=1  k= 1 

~N 

E(4/3)^f 

1=1 

where  Sk  capl  +  Sk  cap2  is  the  total  surface  of  the  two  caps  at  the  inter¬ 
section  of  the  two  particles  and  SktCyi  is  the  surface  of  the  bonding 
between  the  two  particles  (approximated  by  a  cylinder).  Note  that 
the  sum  is  made  on  all  particle  free  surfaces  for  a  given  electrode, 
irrespective  of  their  nature  (NiO,  Ni,  YSZ,  or  LSM). 

Fig.  6  shows  the  evolution  of  the  pore  surface  area  per  unit 
volume  as  a  function  of  the  porosity  e  attained  during  sintering 
for  the  NiO/YSZ  electrode  before  reduction,  without  pore  former 
(NiO/YSZ)  and  with  pore  former  (NiO/YSZ-pf).  Quite  predictably, 
the  pore  surface  area  decreases  as  porosity  e  decreases.  More 


cap2 


(6) 


6(1  -e) 


2  R 


(7) 


when  contacts  between  particles  are  small.  Taking  into  account  this 
particle  size  effect,  we  note  that  our  results  agree  with  the  exper¬ 
imental  data  of  Song  et  al.  [7]  on  LSM-YSZ  composites  sintered  at 
three  different  temperatures  which  lead  to  various  final  particle 
sizes.  These  authors,  using  line  intercept  method  to  analyze  digi¬ 
tal  images  of  their  microstructures,  have  measured  free  surface  per 
unit  volume  of  5, 7  and  1 0  |xm2  |xm-3  for  1  pan,  600  nm  and  400  nm 
approximate  particle  size,  respectively  (their  Fig.  3). 

Fig.  7  shows  the  effect  of  reducing  NiO  to  Ni  on  the  free  sur¬ 
face  for  ONi0  =  65%.  In  accordance  with  the  porosity  increase  due  to 
reduction,  Fig.  7  indicates  that  pore  surface  increases  as  Ni  reduc¬ 
tion  progresses.  Inspection  of  Eq.  (6)  shows  that  the  total  pore 
surface  is  the  sum  of  a  term  that  will  decrease  with  reduction  (sur¬ 
face  of  Ni  particles)  and  of  a  term  which  will  increase  as  contacts 
between  Ni  and  between  Ni  and  YSZ  particles  disappear  (the  cap 
surface  of  the  contacts).  The  third  term  in  Eq.  (6)  can  be  neglected  in 
most  conditions.  Thus,  it  is  not  surprising  that  the  overall  increase 
in  free  surface  is  limited  during  reduction  and  that  its  evolution 
is  not  monotonic.  For  ONi0  =  35%,  we  observed  similar  results  to 
those  shown  in  Fig.  7,  indicating  that  YSZ  percolation  does  not  play 
a  dominant  role  in  the  evolution  of  S. 

Pore  surfaces  were  measured  experimentally  in  the  literature 
for  Ni-YSZ  cermets  by  scanning  electron  microscope  [44,49]  and 
nondestructive  tomographic  imaging  technique  [50].  The  total 
free  surface  reported  in  these  works  are  2.4  pan  [xnrr2 [44,49]  and 
6.8  pm  pm-2  [50]  for  particle  sizes  that  are  somewhat  larger  than 
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Fig.  7.  Effect  of  NiO  reduction  on  the  NiO/YSZ  and  NiO/YSZ-pf  electrodes  for 
d^Nio  =  65%. 

those  used  in  our  simulations.  Again,  these  numbers  are  consistent 
with  the  results  shown  in  Fig.  7  and  with  Eq.  (7). 

5.2.  Triple  Phase  Boundary  Length 

The  total  length  of  TPB  per  unit  volume,  LTPB,  has  been  calculated 
similarly  to  the  pore  surface  by  considering  each  contact  between 
LSM  and  YSZ  particles  or  between  Ni  and  YSZ  particles.  For  com¬ 
pleteness,  we  have  also  calculated  a  TPB  length  between  NiO  and 
YSZ  particles  to  study  the  effect  of  reduction  on  the  geometry  of 
bonded  contacts  between  NiO  and  YSZ  particles.  Ltpb  is  calculated 
by  summing  the  perimeters  2nab  for  those  contacts  (Fig.  3). 

Fig.  8  shows  the  evolution  of  LTPB  with  macroscopic  porosity  s  for 
non-reduced  NiO/YSZ  electrodes  and  for  YSZ/LSM  electrodes.  Since 
the  average  number  of  bonds  and  the  average  bond  size  ab  increase 
as  sintering  progresses,  LTPB  increases  with  decreasing  porosity. 
Flowever,  contrarily  to  Fig.  6,  Fig.  8  indicates  that  it  is  advantageous 
to  introduce  pore  formers  to  increase  LTPB.  For  a  total  given  macro¬ 
scopic  porosity  e ,  the  introduction  of  pore  formers  allows  a  60-70% 
increase  in  LjPB.  Although  verified  only  for  the  NiO/YSZ  electrode, 
this  result  is  likely  to  remain  valid  for  any  type  of  porous  electrode. 
This  is  because  for  a  given  macroscopic  porosity,  the  introduction 
of  large  pore  formers  imposes  a  denser  packing  of  the  remaining 
particles,  thus  increasing  the  average  bond  size  ab  and  the  resulting 
ITPB  [30]. 

Caution  should  be  taken  in  ascribing  too  much  a  positive  role 
to  pore  formers  on  the  increase  of  performance  of  the  porous 


Fig.  8.  Evolution  of  IT pb  with  porosity  for  non-reduced  NiO/YSZ  electrodes  and  for 
the  YSZ  +  LSM  layer  of  the  LSM/YSZ-pf  electrode. 


Fig.  9.  Evolution  of  LTpb  upon  reduction  for  NiO/YSZ  without  pore  former  and  initial 
porosity  £,•  =  0.27  for  two  volume  fractions  of  NiO. 


electrode.  As  shown  in  Fig.  1,  there  are  two  populations  of  pores  in 
microstructures  with  pore  formers:  the  natural  small  pores  left  in 
between  particles,  and  the  large  pores  that  are  left  after  the  removal 
of  fugitive  pore  formers.  The  zones  in  between  large  pores,  which 
contain  many  TPB  points  brought  by  enhanced  packing,  must  still 
be  able  to  diffuse  gases  in  order  to  obtain  good  overall  electrochem¬ 
ical  performances  and  thus  cannot  be  too  dense. 

The  lowest  values  of  LjPB  are  found  in  Fig.  8  for  the  YSZ  +  LSM 
layer.  This  is  due  to  the  already  mentioned  fact  that  LSM  particles 
are  larger  (1  pan)  than  NiO  and  YSZ  particles  (500  nm).  As  for  the 
free  surface  calculated  in  the  preceding  section,  the  calculated  ITPB 
should  be  considered  as  a  lower  bond  since  we  consider  spheri¬ 
cal  particles  as  opposed  to  more  realistic  non-spherical  particles. 
Conversely,  polydispersity  (as  opposed  to  our  narrow  size  distribu¬ 
tion)  should  decrease  the  value  of  ITPB  as  shown  by  Kenney  et  al. 
[21]  who  have  calculated  that  realistic  size  distribution  may  lead 
to  a  decrease  of  LtPB  by  as  much  as  40%. 

The  effect  of  reduction  of  NiO  on  LjPB  is  shown  in  Fig.  9  for  the 
case  of  an  electrode  without  pore  former  and  an  initial  porosity 
Si  =  0.27.  When  the  volume  fraction  of  NiO  is  small  (ONi0  =  0.35),  a 
rigid  skeleton  of  YSZ  particles  is  formed  as  described  in  Section  4. 
The  consequence  is  that  volume  reduction  of  NiO  to  Ni  is  accom¬ 
modated  predominantly  by  pore  opening  in  the  YSZ  skeleton  and 
by  contact  losses  without  macroscopic  shrinkage.  The  consequence 
is  a  decrease  of  the  overall  ITPB  as  shown  in  Fig.  9.  On  the  con¬ 
trary,  when  no  rigid  YSZ  skeleton  is  formed  (ONi0  =  0.65),  volume 
reduction  is  accommodated  by  macroscopic  shrinkage  and  we  have 
actually  observed  that  in  that  case,  new  contacts  between  Ni  and 
YSZ  particles  may  actually  be  brought  by  reduction.  In  that  case, 
the  decrease  in  ITPB  is  much  less. 

All  electrodes,  with  or  without  pore  formers,  exhibit  the  same 
general  trend  described  in  Fig.  9.  This  is  shown  in  Fig.  10a  and  b, 
where  the  variations  inITPB  for  various  initial  porosities  is  collected. 
The  figure  indicates  that  the  increase  in  porosity  is  coupled  with  a 
large  decrease  in  LTPB  for  samples  with  a  rigid  YSZ  skeleton.  Also, 
it  may  be  noticed  that  the  gain  in  LTPB  brought  by  the  use  of  pore 
formers  is  preserved  upon  reduction  for  a  given  relative  density. 
Thus,  it  may  be  concluded  that  pore  formers  should  have  beneficial 
effect  on  ITPB  even  when  accounting  for  the  reduction  stage. 

When  no  pore  former  is  present  and  all  particles  are  of  the  same 
size  2R,  we  have  proposed  a  simple  approximation  for  ITPB  as  a 
function  of  the  composition  of  the  composite  in  one  of  the  particle 
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Fig.  10.  Evolution  of  LTPB  with  porosity  upon  reduction  for  NiO/YSZ  without  pore  former  (a)  and  with  pore  former  (b).  The  effect  of  the  volume  fraction  of  NiO,  ONi0,  is 
demonstrated  for  both  types  of  electrodes. 


type,  (I*  [18]: 

i-TPB  -  3(1  ~2g)Z<Kl  ~  )V3>  (8) 

where  Z  is  the  average  contact  number  per  particle  and  £0  is  the 
initial  porosity  of  the  particle  packing,  before  sintering.  The  exact 
value  of  £0  is  difficult  to  ascertain.  However,  for  the  isotropic  sin¬ 
tering  of  a  layer  which  shrinks  from  an  initial  thickness  to  to  a  final 

thickness  t,  the  term  ((1  -  £0)/(l  -  £))  maybe  estimated  as  t/t0. 
Approximating  Z  by  6,  a  simpler  and  more  useful  version  of  Eq.  (8) 
can  thus  be  written: 

Lm=  ?2(^~e)<l>(l—  <D)^1  —  ~t  (9) 

where  d  is  the  particle  size.  Eq.  (9)  only  introduces  material  and 
process  parameters  that  should  be  accessible  to  experimentalists. 

Eq.  (9)  may  be  compared  to  the  Ltpb  value  of  the  monomodal 
NiO/YSZ  powder  sintered  numerically.  Fig.  8  shows  that,  although 
it  overestimates  LTPB  as  compared  to  the  more  realistic  numerical 
simulations,  it  is  still  a  satisfactory  approximation.  The  overestima¬ 
tion  comes  from  the  approximation,  used  in  deriving  Eq.  (9),  that 
particles  do  not  rearrange  during  sintering  (also  called  the  affine 
assumption).  This  leads  to  an  overestimation  of  ab  and  thus  to  an 
overestimation  of  Ltpb- 

Eq.  (9)  shows  that  decreasing  the  particle  size  is  the  most  effi¬ 
cient  way  to  increase  TPB  length  per  unit  volume.  In  reality,  this  may 
be  achieved,  for  a  given  initial  powder  size,  by  decreasing  the  sinter¬ 
ing  temperature  to  avoid  coarsening  [4].  Furthermore,  for  a  given 
particle  size,  increasing  contact  number  and  decreasing  porosity 
also  allows  larger  LTPB  to  be  obtained.  The  0(1-0)  term  in  Eq. 
(9)  indicates  that  having  a  mixture  close  to  50-50%  composition 
of  ionic  and  electronic  conducting  particles  is  beneficial  for  LTPB 
(although  this  does  not  ensure  good  overall  conductive  properties 
due  to  percolation  effects  as  demonstrated  in  [22]). 

Another  useful  feature  of  Eq.  (9)  is  that  it  explains  quite  simply 
the  large  scatter  of  Ltpb  values  found  in  the  literature.  For  Ni-YSZ 
anodes  quantified  by  coupling  focused  ion  beam  and  scanning  elec¬ 
tron  microscopy  (FIB-SEM),  LTpb  values  have  been  reported  by  three 
different  teams.  Iwai  et  al.  found  LTpb  =  2.6  |xm  |xm-3  at  50%  poros¬ 
ity  and  for  ~3  p,m  (Ni)  and  ~2  p,m  (YSZ)  particles  [51  ].  Wilson  et  al. 
reported  LTpb  =  4.3  [xm  purr3  for  ~2  |xm  Ni  and  ~1  p,m  YSZ  parti¬ 
cles  and  20%  porosity  [44].  Shearing  et  al.  obtained  1 0  |xm  parr 3 
to  13fxm|jim_3  for  approximately  1.5  pun  Ni  and  YSZ  particles 
and  10%  porosity  [52,53].  The  above  reported  values  of  LTpb  are 
consistent  with  each  other  when  considering  the  approximate 


proportionality  rule  derived  from  Eq.  (9): 

f-TPB  OC  L— _ (10) 

Thus,  we  conclude  that  the  large  scatter  found  in  the  literature  on 
Ltpb  values  is  mainly  due  to  porosity  and  particle  size  differences 
between  the  samples. 

6.  Conclusions 

Numerical  microstructures  representing  partially  sintered 
porous  electrodes  have  been  generated  using  discrete  simulations. 
The  originality  of  the  present  approach  is  to  include  the  sintering 
process  in  the  simulation  to  allow  for  realistic  microstructures  to 
be  generated.  We  believe,  in  particular,  that  the  bilayered  struc¬ 
tures  and  the  structures  incorporating  pore  formers  that  have 
been  numerically  generated  here  and  which  are  now  routinely 
encountered  in  real  electrodes,  offer  an  interesting  tool  for  studying 
important  microstructural  characteristics. 

It  has  been  shown  that  the  introduction  of  pore  formers  has 
mainly  a  positive  effect  on  the  microstructure  of  the  electrode.  This 
is  because,  although  it  reduces  slightly  the  surface  area  in  the  elec¬ 
trode,  large  pores  brought  by  pore  formers  provide  an  easy  path  for 
species  diffusion  and  increase  the  TPB  length  in  the  electrode. 

The  present  results  also  show  the  importance  of  the  NiO  vol¬ 
ume  fraction  upon  reduction  of  NiO-YSZ  composites.  When  the 
YSZ  particles  are  forming  a  rigid  skeleton,  reduction  leads  to  a  large 
increase  in  porosity  and  a  decrease  in  the  TPB  length.  When  the  NiO 
volume  fraction  is  large  enough  to  prevent  the  formation  of  a  YSZ 
skeleton,  macroscopic  porosity  does  not  increase  as  much.  Instead, 
macroscopic  shrinkage  is  observed. 

Some  oversimplifications  in  our  model  and  in  the  numeri¬ 
cal  microstructures  should  be  addressed.  Although,  our  sintering 
model  correctly  takes  particle  size  into  account,  it  assumes  that  all 
particles  of  a  given  size,  irrespective  of  their  material  (NiO,  YSZ, 
LSM)  sinter  with  the  same  kinetics.  This  simplification  could  be 
addressed  by  empirically  comparing  real  electrode  microstructures 
with  those  generated  numerically  and  modifying  sintering  param¬ 
eters  (typically  the  surface  energy,  ys)  for  each  type  of  contact  to 
approach  the  real  microstructure.  Also,  the  narrow  particle  size  dis¬ 
tribution  assumed  here  is  clearly  an  oversimplification  that  should 
be  addressed.  The  discrete  element  method  is  a  very  effective  tool 
for  simulating  size  distribution  effects  in  powders  [55,54].  This  is 
certainly  an  avenue  for  simulating  even  more  realistic  porous  elec¬ 
trodes. 

Finally,  the  present  numerical  microstructures  may  be  further 
used  for  mechanical  treatment  to  obtain  the  elastic  and  fracture 
properties  of  porous  electrodes.  In  particular,  they  will  be  useful 
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to  quantify  the  effect  of  large  pores  or  to  evaluate  the  effect  of  NiO 
reduction  on  the  mechanical  properties  of  the  electrodes.  This  work 
is  presently  in  progress. 
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